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ABSTRACT 


This thesis investigates the combined use of ocean models, such as idealized surface current 
flows, and search models, including expanding area and discrete myopic search methods, to 
improve the probability of detecting a near-surface, drifting object over time. Enhanced search 
effectiveness is facilitated by the use of robotic search agents, such as a tactical unmanned un¬ 
derwater vehicle (UUV) or unmanned aerial vehicle (UAV), leveraging simulation methods to 
inform the search process. The presented work investigates the impact of using naive versus op¬ 
timized search patterns on localizing a drifting object, including a surrogate ocean model using 
idealized flow as well as historical data sets with Weibull-distributed perturbations. Numerical 
studies and extensive analysis using different permutations of model parameters (including the 
relative speed of the drifting object, time late in the searcher’s arrival to the search area, sen¬ 
sor sweep width, and duration of the search mission) identify the significant factors affecting 
the overall probability of detection. Such insights enable further explorations using empirical 
datasets for specific oceanographic regions of interest. 
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CHAPTER 1: 

Introduction 


The prevalence of unmanned systems in their diverse application areas is generating not only 
new technologies and capabilities, but is also facilitating new application domains that lie at 
the interface of previously disparate communities. One such interdisciplinary arena is between 
ocean sciences, search missions modeling, and effective employment of autonomous search 
agents such as unmanned underwater vehicles (UUVs) and unmanned aerial vehicles (UAVs). 

In this context, the search for drifting objects in the ocean represents an exemplar of a mis¬ 
sion where optimal execution of the search effort benefits greatly by incorporating information 
about the ocean surface current and environment. Though the respective communities have 
recognized the operational value in this integration, further advances in computation, modeling 
and simulation, as well as robotic sensors are rapidly facilitating its realization. Drifting ob¬ 
jects whose motions are induced by surface currents are relevant to both civilian and defense 
applications, including search for crashed aircraft wreckage (as in Air France 447 [1]), search 
and rescue operations for disabled vessels [2], or clearance missions for near-surface mines 
deployed in commercial waterways [3]. 

In the challenging case of the mine clearance scenario, even if an initial route has been cleared, 
there is no guarantee that the cleared path will not be “re-disturbed” again by the ocean currents. 
In such cases, as well as in search and rescue operations, the time necessary to find one or more 
surface or near-surface drifting mines is a critical measure of performance of any solution strat¬ 
egy. The advent of unmanned systems has created a whole new spectrum for concepts of opera¬ 
tions that allow for more timely deployment of assets and the ability to cover a wide area search 
at diminished expense (either in manpower, risk, or search effort). Any information—either 
historical or sampled—regarding the ocean surface current in the area of operations can poten¬ 
tially aid greatly in informing the search process and enabling more efficient use of searchers in 
terms of minimizing the time to find the mine(s). 

1.1 Motivation 

The Office of Naval Research (ONR) released a Broad Agency Announcement (BAA) in July 
2011 to search for an innovative technology solution leading to wide-area detection, classifi¬ 
cation, localization and tracking of surface and near-surface drifting/oscillating mines and pre- 


1 




diction of mine drift from direct environmental observations and ocean models to support tac¬ 
tical planning (ONRBAA11-024). The envisaged concept of operations, as stated in the BAA, 
utilizes a Mine Drift Tactical Decision Aid (TDA) to predict drift trajectories using real-time 
environmental observations from ocean model data and other sources such as in-situ environ¬ 
ment information retrieved from unmanned systems. The unmanned systems may be equipped 
with sensors such as Side-Looking Sonars, Airborne Laser Mine Detection System, Synthetic 
Aperture Radar, Electro-Optics and so on, to perform wide-area search as well as real-time de¬ 
tection, classification and localization of surface and near-surface drifting or oscillating mines. 
The TDA will also be used to plan the efficacious deployment of the unmanned systems. 

The BAA thus presented an excellent opportunity where my research could complement prior 
work on this topic and provide some insights on the differences in detection probabilities when 
either naive or optimized search pattern is employed in localizing drifting objects in the Persian 
Gulf. The Persian Gulf region was chosen due to its strategic significance in the global oil trade 
as well as the availability of past historic ocean current data that are readily accessible to test 
the hypothesis. In addition, part of the research work has been submitted to the 2012 IEEE/RSJ 
International Conference on Intelligent Robots and Systems (IROS 2012). 

1.2 Economic Significance of Persian Gulf 

Given that petroleum and petroleum-related products remain central to the stability and growth 
of the world economies, the strategic significance of the Persian Gulf cannot be overlooked. 
Two-thirds of the world’s oil reserves and over 40 percent of the natural gas reserves are concen¬ 
trated in the Persian Gulf, a region that is undergoing profound changes both in rapid economic 
and infrastructural development as well as social and cultural changes. In 2007, 17 million bar¬ 
rels of oil a day passed through the Strait of Hormuz (SOH), which is the equivalent of about 
40 percent of all global oil transportation [4]. 

The U.S. has a vested interest in the stability of the Persian Gulf region and has continually 
forward-deployed its Naval forces there for the past thirty years to “provide persistent maritime 
dominance, power projection and effective crisis response” [5]. However, regional tensions 
continue to rise, which may lead to the destabilization of the region. The future of Iraq, the 
Arab-Israeli conflict, asymmetric threats including terrorism, and the role of Iran may once 
again bring the region into chaos [6], 

The events of the last decade have demonstrated that Iran, the preeminent indigenous power 
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in the Persian Gulf, has become a major and indispensable political player in both Iraq and 
Afghanistan, where the United States is militarily involved. Because of this, Washington can¬ 
not formulate policies to these conflicts without taking into consideration Iranian interest [7], 
Notwithstanding, Iran is defiant in its quest for nuclear power and has repeatedly stated that any 
hostile actions taken by the international community, specifically the U.S. or Israel, may result 
in the closure of the SOH. Iran’s threats to disrupt oil have come from as high up as Supreme 
Leader Ayatollah Ali Khamenei. If any country attacks Iran, “shipment of energy from this 
region will be seriously jeopardized,” Ayatollah Khamenei said in a June 2006 speech. Given 
the volatile environment surrounding the Persian Gulf region, the possible closure of SOH due 
to mines is very real. 

1.3 Contributions 

The main contributions of this thesis include the modeling and simulation of time-varying 
search areas that incorporate ocean models, which can be used to inform and improve the search 
for drifting objects by tactical unmanned systems. The search problem formulation, extensive 
simulation studies, and subsequent statistical analysis provide insights into how to incorporate 
real-time and/or historical oceanographic data to maximize the probability of detection. 

1.4 Thesis Organization 

• Chapter 2 provides a glimpse into the world of mine warfare as well as a brief description 
on the geographic features and environmental factors that surrounds the Persian Gulf. 

• Chapter 3 describes the analytic models and relevant works for both search operations 
and ocean surface currents. 

• Chapter 4 investigates the integration and implementation of the model and discusses the 
results and statistical insights obtained from this study. 

• Chapter 5 discusses the implications of the experimental results as well as recommenda¬ 
tions for future avenues of research. 
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CHAPTER 2: 
Illustrating the Dynamics 


Mines are relatively low-cost and highly effective weapons that are used to impose sea control 
by neutralizing or destroying enemy ships as well as interrupting sea lines of communication by 
preventing the passage through an area of sea. The U.S. Navy uses three categories to classify 
mines according to how they are deployed in the water. The categories are bottom, moored and 
drifting mines. Bottom mines have negative buoyancy, sink to the bottom of the ocean and stay 
there. Moored mines have their explosive charge and firing mechanism housed in a positively 
buoyant casing and tethered to an anchor on the sea bottom. The length of cable can be preset 
to ensure that the mines are deployed at the predetermined depth. Drifting mines are mines that 
float freely at or near the surface of the water. 

Psychologically, just the threat that mines exist in an area can easily deter an enemy from 
deploying his naval forces unless the area has been swept by mine sweepers and declared safe 
for transit. Similarly, shipping companies will not want to risk transporting their goods through 
these dangerous waters, and this may result in delays with the shipping schedules as ships may 
be required to make a detour via an alternate route or to call at an alternate port to unload their 
goods. While international law requires nations to declare when an area is mined, the warnings 
need not have to be specific. For example, during World War II, Britain simply declared that it 
had mined the English Channel, the North Sea and the French coast in order to deter the German 
invasion. 

If a rogue nation or terrorist group was able to obtain or build a series of contact mines that 
are currently considered as design relics from World War I and deployed these mines randomly, 
it would effectively choke an entire body of water as no ships will be able to sail through the 
mined waters without a Mine Countermeasure (MCM) vessel to lead the way. If these mines 
were allowed to drift at or near the surface with prevailing currents in the same body of water, 
then they would present an alarmingly difficult and time consuming problem for MCM search 
analysts to track down and dispose of each and every mine. Even if an initial route has been 
cleared, there is no guarantee that the cleared path will not be “mined” again by the natural 
ocean currents. 
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2.1 Brief History of Mine Warfare 

Mine warfare began in 1776 when David Bushnell invented the “Bushnell’s Keg.” These kegs 
were placed in the Delaware River so that they would float into British ships downriver. These 
early mines were called torpedoes. In fact, when Admiral Farragut said “Damn the torpedoes” 
at the Battle of Mobile Bay, he was actually referring to mines [8]. 

When the American Civil War broke out in 1861, the Confederate forces did not have a com¬ 
parable naval force that could match the Union side. Given their limitations in constructing 
conventional coastal defenses, the Confederates chose defensive mining to defend their ports 
and coastal cities. At the end of the war, a total of forty-three Union ships were struck by the 
Confederate mines, twenty-seven of which sank. 

During World War I, mines again played an important role in the naval battle. The allied forces 
sowed a 250-mile wide minefield between Scotland and Norway to stop the German U-boats. 
This minefield was planted in five months during 1918 by American and British ships and 
contained over 72,000 mines. 

During World War II, influence mines (using a combination of magnetic, pressure and acous¬ 
tic sensors) were developed by the Germans. The small number of influence mines laid by 
Germany in 1939 virtually brought the infested ports to a standstill until sweeping techniques 
could be developed. Mines were used by all sides throughout World War II, and surface ships 
were fitted with paravanes to fend off contact mines. In Japan, Operation Starvation saw U.S. 
Army B-29 aircraft sowing more than 12,000 mines around Japanese shipping routes and harbor 
approaches, sinking 650 ships and choking off all maritime shipping [8], 

On the other hand, the U.S. minesweepers and minelayers were primarily used for defensive 
mining operations, protecting friendly harbors from enemy infiltration. Thousands of these 
mines were laid outside Chesapeake Bay, around Cape Hatteras, and at Key West to protect U.S. 
shipping against enemy submarines. In addition, large mine fields were laid in the Atlantic, off 
the coasts of Trinidad and North Africa as well as in the Mediterranean. While there is no record 
of any enemy ships being sunk or damaged in these mine fields, neither is there any record of 
any enemy ships passing through these fields. 

It was perhaps the knowledge that the mine fields existed in these areas in the first place that 
deterred the enemy from making any attempts to transit through these areas at all. During 
the Korean War, American forces were prevented from landing at Wonsan for over a week in 
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1950 until 3,000 Communist mines could be cleared. Similar tactics were used again during 
the Vietnam War, when the Vietcong planted numerous mines made from modified bombs in 
shallow rivers and deltas to disrupt shipping. 

Throughout the history of the U.S. Navy, there were more ships damaged or lost to mines than 
missiles [9]. During the Korean War, mines laid by North Korean forces damaged 11 U.S. naval 
vessels. During the Iran-Iraq War, an Iranian M-08 mine costing $1500 was able to put a 25- 
foot (8 m) hole in the hull of the USS Samuel B. Roberts (FFG-58), forcing the ship to seek 
temporary repairs in a dry dock in Dubai, UAE and causing a whopping 96 million dollars in 
damage [10]. During the Persian Gulf War, both surface and air mine countermeasures forces 
were deployed to find and neutralize primitive Iranian mines. Two U.S. navy vessels (USS 
Princeton [CG-59] and USS Tripoli [LPH-10]) were severely damaged by Iraqi mines laid in 
areas previously thought to be safe for transit. 


2.2 What are Drifting Mines? 

Mines come in two general varieties: contact and influence. For detonation, contact mines 
require a physical impact between the target and the mine itself, and are either moored in the 
mid-water region to an anchor on the bottom, or simply set adrift until they encounter a target. 
A drifting mine is simply a floating mine without any mooring. These mines usually weigh 
about 120 kg, including 80 kg of explosives (TNT). 

Drifting mines were occasionally used during World War I and World War II. However, they 
were more feared than effective. After World War I, the drifting contact mine was banned, 
but continued to be used during periods of sea conflict. The drifting mines were much harder 
to remove after the war and they caused about as much damage to both sides. According to 
the terms of the Hague Convention, when mines break away from their moorings and become 
drifting mines, their design should deactivate the arming mechanism and render the mine inert. 

However, in practice, one cannot take this for granted as there is possibility that after several 
years at sea, the deactivation mechanism might no longer function as intended. In such situa¬ 
tions, the mines may become active again and pose a danger to merchant shipping as well as 
civilians when swept ashore. Such danger still persists in present time and the UK Department 
for Transport issued an updated notice (Marine Guidance Note MGN 323 (M+F), “Explosives 
Picked Up At Sea”) to all mariners in 2006 on the hazards that these explosives pose. 
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2.3 Offensive Mine Warfare 

Mine warfare can be classified into offensive and defensive mine warfare. Forces that use offen¬ 
sive mine warfare always have the initiative and the upper hand while the mine countermeasures 
groups are the underdogs who are constantly seen as responding to those initiatives. The fol¬ 
lowing paragraphs will take a look at three types of offensive mine warfare that a rogue nation 
or terrorist groups may employ against a superior naval force like the U.S. 

2.3.1 Interdiction 

In mine warfare terms, the use of mines for interdiction involves the obstruction of shipping 
lanes and transit points with the intention to restrict the ability of hostile forces to mount attacks 
from the sea or to inhibit the free flow of commerce or supply routes [11]. In many cases of 
interdiction mining, the minefields are covertly laid off enemy ports or assets. Once such a 
minefield is discovered, the port in question is sealed off until all the mines have been swept. 
It is thus a massive effort required of the mine countermeasures vessels to carry out routine 
surveillance at every port, on all shipping routes, all of the time. 

2.3.2 Attrition 

Attrition can be defined as the art of mining High Value Units into a selected area or prevent¬ 
ing them from getting to sea. Examples of HVUs are an aircraft carrier or a ballistic missile 
submarine. Usually attrition is achieved by selecting a suitable choke point through which the 
High Value Units need to pass and then mining it. In the world today, many nations have a 
submarine force of two to three boats. If these boats are unable to get out of port, then all the 
investment placed in acquiring, supporting and maintaining that force would have been wasted 
and that nation’s submarine force is effectively nullified by mines. 

2.3.3 Terrorism 

A minefield sowed by terrorist groups can consist of as little as a single covertly placed mine or 
random sowing of mines in a trade route. The terrorist’s main objective is not the destruction 
of ships but the economic damage that follows as a result of disruption to trade and the political 
fallout between the government and its citizens for failing to keep the waters safe. The business 
community and the general public, suffering from the economic devastation of the trade dis¬ 
ruption, may then add pressure on the government to change their laws or policies, which will 
eventually favor the terrorists. 
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2.4 Geographic Features of Persian Gulf 

The Persian Gulf, also known as the Arabian Gulf, is a shallow, semi-enclosed basin situated in 
an arid zone where evaporation greatly exceeds precipitation and river runoff. The Persian Gulf 
is approximately 990 km long, with a minimum width of about 56 km and maximum depth of 
100 m near the SOH, a maximum width of about 338km, and a mean depth of about 35m in 
the middle of the gulf. The inland sea occupies an area of approximately 251,000 km 2 . River 
inflows are mostly in the northern end of the gulf, primarily on the Iranian side, with the largest 
being the Shatt-Al-Arab, a river formed by convergence of the Tigris, Euphrates, and Karun 
rivers. 

The gulf is surrounded by the deserts of the Arabian Peninsula to the west and the Zagros 
Mountains of Iran to the east (shown in Fig. 2.1). Iran’s Zagros Mountains begin in northwestern 
Iran and end in the southern province of Hormozgan, north of SOH. These mountains act to 
channel the winds to the southeast. Thus, the Persian Gulf is characterized by northwesterly 
or westerly winds throughout the year with mild and continuous summer winds as well as 
sudden and violent winds in the winter that are associated with synoptic weather systems [12]. 
Furthermore, the shallow water in the Persian Gulf means that the difference between the wind- 
driven current at the surface and sea bottom is minimal and the current flow can be approximated 
by two-dimensional flow plots [13]. This thesis will focus on the northern Persian Gulf region 
that is separated from the southern region by a front off Qatar. 


2.5 Environmental Factors From Previous Studies 

There are three primary contributors to mixing in the water column: (1) tides, (2) winds and 
waves, and (3) evaporation. Internal waves at density interfaces, and topographic features such 
as islands, also contribute to mixing, but their contribution is unknown and probably small [12]. 
The kinetic energy of the water velocity as a result of these factors can be partitioned among the 
three primary contributors in the approximate ratio of 100, 10 and 1, respectively [12]. However, 
each factor has a different scaling time: tides vary every few hours at diurnal periods, wind- 
driven currents vary over a few days, and density-driven currents or thermohaline circulation 
varies over a period of weeks. 

The strong tidal currents create a bottom, mechanical, turbulent friction layer. The enhanced 
turbulent mixing homogenizes the water column from the bottom up, and a bottom mixed layer 
is evident in all hydrographic sections. The other two mechanical mixing processes are strongest 
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Figure 2.1: Topographic Map of the Middle Eastern countries centering on the Persian Gulf. 



at the surface. Wind friction creates a surface shear layer and also creates surface wave fields, 
both of which form a surface mixed layer. 

The currents in the northern portion of the Persian Gulf are dominated by northwesterly winds 
that vary in speed from 10 to 15 knots. Together with freshwater discharge from the Shatt-al- 
Arab, this creates a southeastward-flowing coastal current along the Arabian and Iranian coasts 
as shown in Fig. 2.2. Residual currents indicate a surface inflow of the Indian Ocean Surface 
Water into the Persian Gulf as well as a deep outflow of the Persian Gulf through the SOH. 

During winter (January-March), the surface water flowing along the Iranian coast seems to 
form an anticlockwise gyre or basin-wide cyclonic circulation in the southern half of the Gulf. 
With the onset of the southwest monsoon in summer (July-September), the current becomes 
progressively stronger and penetrates far into the northern gulf because of the strengthening 
stratification, increasing inflow and weakening winds. The evaporation of fresh water enhances 
mixing by increasing the salinity and density of the surface water, thus reducing surface stability. 
The combined cooling and evaporation increases the density sufficiently to overturn the water 
column and creates a well-mixed water column north of Qatar. The circulation in the Persian 
Gulf is therefore highly dependent on the reversal of the monsoon [14]. 

Other factors that may affect water circulation in the Persian Gulf are the Coriolis force, river 
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runoff, bottom topography and coastline boundaries. The effect of the Coriolis force is to deflect 
the surface inflow towards the Iranian coast and the subsurface outflow towards the Arabian 
coast. The influx of fresh water is seasonal and localized close to the water sources at the 
northwest region as well as along the Iranian coast [15]. 



Figure 2.2: Sketch of surface currents and circulation processes in the Persian Gulf and western 
Gulf of Oman (From Reynolds, 1993). The dashed lines in the central part of the gulf indicate the 
approximate seasonal location of the salinity front between the inflowing fresh water and the more 
saline waters of the gulf. 


2.5.1 Tides 

The tides in the Gulf are complex standing waves and the dominant pattern varies from be¬ 
ing primarily semi-diurnal to diurnal. The tidal range is large, with values greater than 1.0 m 
everywhere [12]. Tides are important for stirring and mixing waters vertically and on a hori¬ 
zontal scale of 10 km, but tides are not an important contribution to the overall circulation in the 
Persian Gulf. The tidal influence was found to be minimal (less than 2 cm/s), except in a few 
localized areas (e.g., the Iranian coast and the SOH) where the tidal residuals are large enough 
to enhance the dominant density-driven flow [16], Tidal current (averaging longer than a day) 
has negligible residual energy, and as a result, basin-scale advection from tides is not consid- 
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ered by most circulation models. Tides are, however, important on smaller scales of horizontal 
length (< 10 km) and time (< 24 hrs). 


2.5.2 Winds 

Low-level winds in the Persian Gulf are dominated by land-sea breeze circulation that varies 
seasonally and diurnally. In summer, the sea breeze is deeper and wider than in winter while on 
a diurnal scale the intense temperature difference between the land and water surfaces during 
the daytime and night causes the land-sea breeze circulation to add a landward component to 
all winds. Over the north coast, the opposing ambient wind creates a sea breeze front that has 
very limited landward penetration. At the south coast, the penetration is over 250 km [17]. 

A study by Zhu and Atkinson showed that in January, wind direction at Kuwait International 
Airport (KWI) was predominantly northwesterly. In January, the direction of the mean wind 
was clearly from the northwest, but on some occasions the wind blew from the east. In April, 
the direction of the daytime mean wind was from the east and southeast, i.e., from the sea. In 
July, the mean direction was again from the northwest throughout the day. The October pattern 
was rather similar to that in April. 

An analysis of each month’s characteristics was carried out by Williams using the wind dataset 
obtained from Fleet Numerical Meteorology and Oceanography Center (FNMOC) and the mean 
wind speed of about 5 m/s as shown in Fig. 2.3 was shown to flow along the major axis [18]. As 
the winds progress down the major axis of the Persian Gulf, their magnitudes tend to increase 
due to the decreased friction from the land-sea breeze circulation winds and the downslope 
wind flow off the Zagros Mountains to the east. 

April’s mean wind currents shows the winds flowing predominantly along the major axis with 
the winds traveling along the Zagros Mountains converging and enhancing the flow in the Per¬ 
sian Gulf as shown in Fig. 2.4(a). With the onset of the southwest monsoon, there is a strong 
area of convergence near the SOH with the majority of the flow moving towards the Arabian 
Peninsula. In the northern Persian Gulf, the areas of weak winds also increase to moderate flow 
as illustrated in Fig. 2.4(b). One thing to note is that both figures indicate the presence of strong 
offshore winds. 

2.5.3 Density 

The Persian Gulf Water is one of the most saline water masses in the world. It is formed 
due to excessive net fresh water loss to the atmosphere through evaporation and constricted 
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Timeseries of the Amplitude of Winds 



Figure 2.3: Plot of Wind Amplitudes from Feb to Jul 2006 (From Williams, 2007) 


Mean Wind Vectors for APR 




Mean Wind Vectors for JUL 


(a) (b) 

Figure 2.4: Illustration of the mean wind vectors for the months of Apr and Jul 2006 as obtained 
from FNMOC. The northwesterly winds predominantly flow along the major axis while the onset of 
the southwest monsoon causes a convergence near the SOH (From Williams, 2007) 


water exchange with the open ocean through the SOH due to the much shallower depths of 
the basin, wider strait to the open ocean and lack of a distinct sill as shown in Fig. 2.5. The 
rate of vertical mixing is highest in the SOH where numerous topographic features constrict the 
channel, accelerating tidal flows and generating large eddies capable of scouring the seafloor. 

The dense Persian Gulf Water exits the gulf as a deep outflow and induces an inverse estuarine- 
type water exchange through the SOH. The low salinity Indian Ocean Surface Water spreads 
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Figure 2.5: Physical characteristics of Persian Gulf. (National Geophysical Data Centre, 2010) 

into the interior of the basin generally in two pathways: one branch (northern branch) propagates 
along the Iranian coast and in summer can reach the basin west of Qatar to form a cyclonic 
circulation there. It drives a returning coastal current along the Arabian coast north of Qatar. 
The other branch (southern branch) veers southward to cross the basin to reach the southern 
UAE coast and also forms an overall cyclonic circulation in the southern gulf [19]. 

In summer, the vertical mixing is suppressed by the stable thermocline and penetrates to a depth 
of about 30 m, allowing the Indian Ocean Surface Water to extend farther into the gulf without 
being greatly modified. In winter, the strong vertical mixing induced by the heat loss leads to an 
upward salt flux that causes the water column to be well mixed vertically to a depth of about 70 
m [20], The densest water forms during winter in shallow water at the northern end of the gulf. 
High salinity water forms along the western and southern Arabian coastlines, but the highest 
densities are not observed there because winter temperatures are milder than those in the north 
and considerable dilution occurs before high salinity water formed in shallow bays reaches the 
main gulf basin [21]. 
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CHAPTER 3: 
Model Development 


3.1 Search Models for Drifting Objects 

The mathematics of searching for objects in motion has an extensive theoretical heritage in the 
operations research community [22], and models for addressing expanding search areas as a 
result of moving objects (such as a fleeing submarine) have been well-studied in contexts called 
the “flaming datum search” problem [23]. Additionally, numerous efforts have recently exam¬ 
ined the employment of robotic search agents, such as unmanned aerial vehicles equipped with 
sensors, including works by [24] and [25]. Among others, these works highlight the advantage 
that even a myopic search strategy of greedily maximizing the detection probability can provide 
in improving the time efficiency of the search. This section briefly describes these two search 
methodologies that are relevant to the desired search for drifting objects. 


3.1.1 Expanding Area Search 

In search problems where the search area grows with time, the challenge of conducting a search 
in this expanding area becomes evident as the searcher often cannot “keep up” with the rate of 
area expansion. An example of this is an instance of the flaming datum search problem, where 
an otherwise concealed submarine makes its location known (e.g., by shooting a missile) at 
an instant in time before submerging and fleeing from that location. Given the uncertainty in 
the submarine’s heading in the 2D plane, coupled with the time necessary to get a search asset 
on station (known as the time late and denoted t), the area where the submarine could possi¬ 
bly be is encapsulated by a. farthest-on circle. This area grows (assuming constant maximum 
speed, U, of the target) quadratically in time, whereas the searcher’s ability to sweep this area 
is only linear. As such, the probability of detection, Pd, of finding the target as a function of 
time can be approximately bounded between pf in(lom < P D (t) < P^ eal , where F™dom anc j pKleal 
are the cumulative detection probabilities when using a random and ideal coverage strategy, 
respectively [26]: 
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where V is the searcher speed and W is the sensor sweep width. 

Examples of ideal search patterns are lawnmower or spiral patterns that methodically cover the 
search area. Although these analytic models are useful for gaining an intuition of the relevant 
search parameters, in realistic conditions where the expanding area is not a symmetric shape, 
e.g., due to ocean currents and eddies, other approaches must be utilized to compute the time- 
varying area expansion. 

3.1.2 Discrete Myopic Search 

The aforementioned random and ideal search strategies do not use prior information about the 
spatial likelihood of target presence to inform the search path. In contrast, methods such as the 
discrete myopic search approach [26] leverage such information to good advantage. Discrete 
myopic search is used to allocate opportunities to conduct a “unit of search” or a look. The 
associated algorithm maximizes the probability of detecting a target by always looking in the 
region that yields the highest increment in probability of detection. For discrete looks, a locally 
optimal plan that maximizes the probability of detecting the target can be shown to result in 
a locally optimal plan that minimizes the mean cost (e.g., time) to find the target [27], which 
makes it a viable approach in the search for time-sensitive objects. 

A given search area (though possibly expanding) is first discretized into n cells in which a 
probability that a target is within the i th cell, denoted p,. is assigned. A budget of K looks is 
available for allocation among these cells. Each successive look in cell i has a probability /3, of 
not detecting the target, with each look assumed independent of all other looks. Given k[ looks, 
for ki G {0,1,... ,K}, in the i th cell, the objective is to maximize the probability of detection as 
defined by: 

= 1-#'■), (3-3) 

i=i 

subject to nonnegative integer decision variables, ki, and a constraint on the total number of 
looks, i.e., Yli= i ki < K and £" =1 pi = 1 [26]. 

The myopic search algorithm outlines a procedure for allocating each look, for which the k+ 1 st 
look, is assigned to the cell with the highest increase in the overall probability of detection, and 
as such, is a “greedy” algorithm (i.e., steepest gradient ascent). Global optimality is ensured, in 
the absence of false positive detections, due to the concavity of the objective function in ki [26], 
We leave for future study the impact of false alarms on the searcher’s effectiveness [28] (e.g., if 
the searcher is penalized for not classifying the mine-like object) and the use of more advanced 
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search strategies that take into consideration the inherenet uncertainty in sensor performance 
(e.g., multiscale adaptive search algorithm for single or multiple searchers [29]). 

3.2 Ocean Surface Current Models 

3.2.1 Computational Ocean Nowcast/Forecast Models 

Current computational ocean models generally use a number of methods (e.g., climatology, 
satellite imagery, and in-situ measurements using sensor floats and buoys) to improve their ac¬ 
curacy in predicting the ocean environment, including its temperature, salinity, and, important 
for our study, surface current speeds. For example, combining satellite data assimilation with 
existing ocean nowcast/forecast systems has shown to provide a better estimate of sound veloc¬ 
ity profiles, which depend on the temperature, salinity, and pressure of the water column [30], 
The Persian Gulf is the region of interest to both science and military communities and has been 
heavily studied and modeled [12,31]. 

Accessible from the National Oceanic and Atmospheric Administration’s (NOAA) website, the 
Navy Coastal Ocean Model (NCOM) is the Navy’s real-time operational global nowcast/forecast 
system developed by the Naval Research Laboratory that is built on the base of the Princeton 
Ocean Model (POM) [32]. The NCOM produces forecasts at three-hour intervals out to 72 
hours in the future. Atmospheric forcing is provided by the Navy Operational Global Atmo¬ 
spheric Prediction System (NOGAPS) while the Ocean Surface Current Analyses - Real time 
(OSCAR) model provides near-real-time global ocean surface currents derived from satellite 
altimeter and scatterometer data. These models provide a data resolution ranging from 1 de¬ 
gree to as high as 1/3 degree increments. An example of a surface current plot generated from 
OSCAR is shown in Fig. 3.1. 

3.2.2 Simplified Ocean Flow Modeling 

The extreme complexity of the computational models described above arises from the immense 
number of state variables and parameters necessary for high-fidelity models. However, the long 
time necessary to generate surface current predictions using these models limits their value in 
tactical or short-time horizon scenarios. Specifically, in the context where we seek to minimize 
the time late to execute the search, time delay for these forecast models may render the search 
for drifting objects untenable. 

Alternatively, we consider the use of analytic models from two-dimensional potential flow anal¬ 
ysis (e.g., for incompressible flow in fluid dynamics) that provide closed form expressions for 
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Mean -> 1 .0 meter/sec (0.514 m/s = 1 knot) 

Figure 3.1: Five-day mean of the global ocean surface currents measured in meters per second. The 
plot is produced by OSCAR and the data is centered on December 26, 2011. The plot covers the 
world oceans from 60°S to 60°N with a resolution of 1°. 

various elementary flow types, such as uniform flow, sources, sinks, vortices, and their various 
superpositions. These potential flow models can be used and solved in computationally expe¬ 
dient manners (i.e., by solving Laplace’s equation) to provide a surrogate for sophisticated but 
expensive computational predictions. Additionally, the flexibility of potential functions is such 
that they can readily be used to construct an approximation to empirical historical data, which 
can then inform the search process in the methods presented in this thesis. 

As an illustrative example, consider a simulated surface current generated by a uniform flow in 
the presence of a generalized source-sink pair illustrated in Fig. 3.2. The ideal flow field lines 
identify the paths (assuming no noise) that a drifting particle would take for arbitrary initial 
conditions. For example, for an object dropped at location marked “1,” its subsequent trajectory 
will follow the darkened path in time. 

3.2.3 Particle-based Expanding Area Simulation 

The above potential flow models serve as computationally efficient surrogate surface current 
models, and can be used to simulate trajectories of a cloud of particles moving in time through 
the field. However, the idealized flow models do not capture the stochasticity present in real 
ocean currents. To address this determinism, we can introduce perturbations of the particle 
paths at each time step, such that the noise causes random dispersions of the particle cloud. In 
this manner, Monte Carlo simulations using a large number of particles can be used to rapidly 
approximate the expanding area in which a drifting object could be present, by taking the convex 
hull of the point cloud as it evolves in time. 
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Figure 3.2: Illustration of simulated uniform potential flow at 1.0 m/s with a generalized source-sink 
pair, which acts as a bluff body in the flow. An object is dropped at (x,y) coordinate (0.2,1.0) and 
will follow the u — v velocities (i.e., in flow relative coordinates) created by the composite flow, shown 
by the magenta line over an evolution of five minutes. 


In the simulation model presented in this thesis, 500 particles are initialized in the flow field. 
As time evolves, Weibull-distributed perturbations are added to the particles’ positions as their 
trajectories are simulated using the potential flow velocities. The primary motivation for the 
use of the Weibull distribution is that it has been shown [33] that the probability distribution 
function (PDF) of the surface current speeds generally follows a Weibull distribution. Based 
on the simulated dimensions, the selected Weibull distribution parameters include a scale factor 
A = 1.0 and a shape parameter, y = 1.5 for each coordinate (See Fig. 3.3), which appear to 
produce noise distributions consistent with the previous study. 

Fig. 3.4 depicts the resulting particle cloud evolutions in the absence of translational flow for 
different configurations of the Weibull distribution parameters. These differences cause biases 
in x and y dimensions, leading to different spatial dispersions. Note that the position of every 
particle is translated by the centroid coordinates of the point cloud to represent simulated motion 
in all coordinate (including negative) directions. 

The enclosing convex hull is formed by the (red) lines, which corresponds to the expanding 
area, which can then be used in place of the circular expanding area formulas described in the 
previous section. 
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(a) (b) 

Figure 3.3: Weibull( 2 -Parameter) Plots for different scale factor A and shape parameter 7 (or referred 
to as k in the plots) parameters (a) Probability Distribution Function and (b) Cumulative Distribution 
Function (Wikipedia, 2012 ) 


3.3 Integrated Search and Ocean Modeling 

3.3.1 Simulation Scenario Setup 

Assume the region of interest where search for drifting objects is in relatively shallow water, 
such that the (most likely wind-driven) current flow is adequately approximated by 2D flows 
at the surface. Consider a 200 km by 200 km region of interest where search operations are 
likely to occur. The scenario commences when a report is received identifying the known time 
and initial location of a drifting object (e.g., last known position of a downed aircraft, human 
intelligence of a near-surface mine deployment). At the time of receipt of the report (which may 
already reflect delayed information), a tactical unmanned aerial sensor platform is deployed to 
conduct search, but has to contend with additional time late before arriving on station. The 
searcher will employ its given search strategy, and begin execution of the search upon arrival to 
its designated start point. 

The objective of the investigation presented in this thesis can now be stated as: 


Given a (surrogate) surface current flow model and its resulting prediction of the 
expanded search area for a drifting object with known initial position, how does 
the choice of search strategy affect the probability of detecting the object? 
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(a) 


(b) 


Figure 3.4: Illustration of the simulated dispersion of particles in the presence of two-dimensional 
perturbations over time. These noises are drawn according to a Weibull distribution, which has 
been shown to model surface current randomness. Depending on the distribution parameters, one 
can construct (a) a balanced, nearly isotropic dispersion pattern, or (b) a biased dispersion along a 
preferred dimension, to best model the ocean environment in the region of interest. 


3.3.2 Search Strategy Comparison 

The two search strategies to be investigated can be characterized as follows. 

Linear Search 

The searcher transits directly to the predicted location (based on time late, t) of the drifting 
object under the influence of idealized flow with no random perturbations, and proceeds to 
execute a circular spiral out pattern for the duration, T, of the search mission. This naive 
approach reflects current standard operating procedures for responding to a search cue, and as 
such, does not account for the ocean model-dependent expanding area nor the particle-based 
distribution information. 

Myopic Search 

The searcher uses the flow model with noise presented previously to predict the course of 500 
simulated particles, recording the convex hull and centroid location information at each time 
step up through t. Once the convex hull has been appropriately computed (which has fHn logo) 
complexity in number of particles n ), the searcher divides this area into four partitions along the 
direction of flow as well as in the orthogonal axis (as in Fig. 3.5). The number of particles in 
each of the four partitions is used to determine a prior probability of the object’s likelihood to 
be in that partition. The searcher then can use the discrete myopic search algorithm described 
in Section 3.1.2 to allocate search effort in one-hour units up to search duration, T, among the 
four partitions. (We can assume the searcher employs an ideal search once in a partition.) 
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Figure 3.5: Simulation of a particle cloud is used to define a time-varying search area (i.e., the convex 
hull enclosing the particle cloud) as it evolves in time. The blue dots represent the locations of the 
500 particles at a time late of 36 hours, and the red dotted outline shows the growing search area 
over each time interval. This area of uncertainty is then decomposed into four partitions, for each of 
which a prior probability of object presence is computed using the percentage of the total particles 
contained within. 

Table 3.1: List of five factors used in the simulation runs to investigate their impact on search 
effectiveness 


Factors 

Typical Values 

Description 

Searcher speed ( V ) 
MQ8 Fire Scout 

50.0 m/s 
« 100 kts 

Search speed of UAV 

Drifter speed (U) 

0.1 to 1.0 m/s 

Min and max sea surface 
current speeds 

Time Late (t) 

12 to 36 hrs 

Time delay before 

Searcher can be on station 

Search Time (T) 

1 to 4 hrs 

Amount of time Searcher 
has in area of interest 

Sweep width (VF) 

100 to 1000 m 

Nominal sensor performance 


3.3.3 Relevant Factors and Design of Simulation Experiments 

Though there may be numerous variables relevant to ocean and search models, the work pre¬ 
sented in this thesis highlights a set of factors believed to be of operational concern for the 
search mission. These factors and their notional values are provided in Table 3.1. 

We treat the initial location of the drifting object as an uncontrollable factor [34], and as such, 
the simulation code (written in Matlab) has been designed to enable such investigations of the 
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effect of arbitrary starting locations, as well as to incorporate other elements such as historical 
surface current data. 

Given that each factor has a range of possible values, exhaustive examination of the factor space 
is not feasible, especially as the Monte Carlo nature of the particle cloud evolution imposes 
additional computational overhead. Overcoming this challenge is the essence of the field of 
statistical experimental design, where well-chosen combinations of factor values can be used 
to statistically infer the relationship between them and the response variable or the measure 
of performance [34]. For example, screening designs of experiments can use permutations of 
upper, lower and mean factor values to extract which factors, e.g., a subset of those listed in 
Table 3.1, are most significant in affecting the probability of detection. Numerical studies and 
their results for the simplified potential flow model highlighted in Section 3.2.2 as well as the 
studies and results based on the historical surface currents data obtained from the Shallow Water 
Assimilation Forecast System (SWAFS) are highlighted in the next chapter. 
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CHAPTER 4: 
Analysis of Model Output 


In designing a simulation model, there is a need to ensure that the various building blocks are 
working in the right manner before we proceed to expand the analysis to include historic data 
sets, where a number of factors may make it difficult to troubleshoot where the model has gone 
wrong. The approach used in this thesis was to perform the search analysis with a simple 
potential flow model first, and then expand the analysis to include historic data sets obtained 
from SWAFS. 


4.1 Search Analysis Using Potential Flow Model 

During the search process, if the drifting object falls within the range of the searcher’s sweep 
width, it is defined as a successful detection and the search is terminated and a value of one is 
assigned to the simulation run. If the searcher is unable to locate the object within the given 
search time, then a value of zero is assigned to the simulation run. The percentage of successful 
detections is obtained by counting the number of runs where the outcome of the simulation is 
set to one over the total number of simulation replications. 

As described in Section 3.3, we chose a number of values from each factor for our simulation 
runs so as to infer the relationship between the search performance and the response variable. 
In addition, we note that the drifter speed is always proportional to the searcher speed and a 
better representation would be to indicate the drifter speed as a ratio of the searcher speed. That 
said, we fixed the searcher speed at 50 m/s for all our analysis and if we assume that the drifter 
is traveling at a constant speed of 1.0 m/s, we explain the analysis by saying that the drifter 
is moving at U = 2% of the searcher speed. Then, we chose the values listed in Table 4.1 
and perform the simulation runs on 500 samples for both the linear search and myopic search 
strategies. Next, we repeat the entire simulation for a different start location of the drifter 
and verify if both sets of results are generally consistent. The total number of simulation runs 
completed for the potential flow model amounts to 144,000 simulation runs. 

Though the data mining and analysis portion highlights various insights, we present the results 
for four analyses where the impact of the search and ocean model parameters on the probability 
of detection objective is evaluated. In the charts below, the horizontal axis denotes maximum 
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Table 4.1: List of values used in the simulation runs to investigate their impact on search effectiveness 
in the potential flow model for both myopic and linear search strategies 


Factors 

Values Used in Simulation 

Drifter speed (U) 

0.2%, 1% and 2% of V; (0.1, 0.5 and 1.0 m/s) 

Time Late (t) 

24 and 36 hrs 

Search Time (T) 

1, 2, 3 and 4 hrs 

Sweep width (W) 

100, 500 and 1,000 m 


duration of the search mission, T, with the dependent axis representing the probability of de¬ 
tection (computed as the percentage of successful simulation trials). The search performance is 
shown for both search strategy types. 

4.1.1 Analysis 1: Time Late, t = 24 hrs; Sweep Width, W = 500 m 

Based on a time late of t = 24 hours, and a searcher with a sensor sweep width W = 500 meters, 
the simulation results for the case where the ratio of the object speed (based on current flow) to 
the searcher speed is U = 2% is shown in Fig. 4.1(a). The linear search approach continues to 
yield improved performance as the amount of time spent searching increases, whereas myopic 
search appears to plateau after two hours of search effort. This is because the drifting object 
has already moved out of the search zone within those two hours such that continued searching 
in subsequent areas is fruitless. We observe a crossover point where, although the myopic 
approach initially excels, linear search begins to perform better after a search duration time of 
about two hours. 

For U at 1% (shown in Fig 4.1(b)) or less, both myopic and linear search have similar per¬ 
formance results that show detections occur nearly 100% of the trials. This fact shows that 
for small enough relative flow speeds, the expanding area does not increase too rapidly to be 
completely covered by a searcher with a sweep width of at least 500m. 

4.1.2 Analysis 2: Time Late, t = 24 hrs; Sweep Width, W = 100 m 

When the sensor sweep width is reduced from 500m to 100m, there is a sharp drop of about 
60% in the detection probabilities for U = 2%, as illustrated in Fig. 4.2(a). The same crossover 
is seen at search times near two hours, where linear search begins to outperform the myopic 
search with increasing T. On the other hand, once again, the performance of myopic search 
tapers off at about T = 2hrs. 
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(a) 



Figure 4.1: Probability of detections from simulation studies comparing linear search (red squares) 
and myopic search (blue diamonds). The search and ocean parameters are t = 24hrs, W = 500m, 
and object speed ratio is set to (a) U = 2%, which sees linear search improving with increased search 
times, and (b) U = 1%, where either search method yields nearly perfect detections due to the short 
time late and large sensor sweep width. 


An interesting outcome is observed when U is reduced to 0.2%, shown in Fig. 4.2(b). There is 
a reversal in performance such that myopic search performs better when more than two hours 
are available for the search mission duration. For instance, myopic search achieves nearly 90% 
successful object detections at T = 4hrs while linear search only manages an object detection 
rate of about 70% for the same duration. Further simulation studies for the range of U values 
show that this transition (between Fig. 4.2(a) and (b)) occurs when U is less than approximately 
1%. Conversely, for faster surface currents, i.e., U approximately greater than 1%, linear search 
serves as the better search alternative, all other factors remaining constant. 



Figure 4.2: Search mission results for t = 24 hours, W = 100 m, and for (a) U = 2 %, and (b) 
U = 0.2%. For shorter sensing distances and slower surface currents as in (b), the myopic search 
approach appears to vastly outperform the naive linear search. 
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4.1.3 Analysis 3: Time Late, t = 36 hrs; Sweep Width, W = 500 m 

Next, we investigate the impact of increased time late of T = 36 hours (with W = 500m), which 
nominally only results in a larger search area due to further expansion over time as compared 
with the above analyses. In doing so, for surface current speed ratios U = 2% (depicted in 
Fig. 4.3(a)), linear search has about a 20% advantage in detecting the drifting object over myopic 
search as early as the first hour of search. 

However, this advantage is quickly lost, as myopic search is able to produce similar search 
performances when search time allocated to the mission is increased to two hours or more as 
illustrated in Fig. 4.3(a). When U is reduced to 1%, as shown in Fig. 4.3(b), the most noticeable 
change is the diminished gap in performance (i.e., to only ~ 10%) between the two search 
strategies at the one hour search mission mark. Note that thereafter, myopic search has a slightly 
better performance as compared to linear search. 



Figure 4.3: Search mission results for t = 36 hours, W = 500 m, and for (a) U = 2%, and (b) U = 1%. 
For faster surface speeds, the myopic search approach suffers initially by too narrowly focusing its 
search effort, although for both speed ratios, it appears to close the gap in search performance for 
missions extending past the initial first hour. 


4.1.4 Analysis 4: Time Late, t = 36 hrs; Sweep Width, W = 100 m 

As before, we examine what happens when the sensor sweep width is reduced to W = 100m, 
for extended time late T = 36hrs. As this is the most challenging operational case, we anticipate 
diminished search performance relative to the previous studies. As shown in Fig. 4.4(a) and (b), 
for U = 2% and U = 0.2%, respectively, we observe that the highest probability of detection for 
either linear or myopic search strategies in the former analysis does not reach 30%, even with 
an extended search mission duration of four hours. This trend is similar to the chart shown in 
Fig. 4.2(a), which similarly had a short sensor range trying to search for faster-moving objects. 
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For slower surface currents, with U = 0.2%, comparing the linear search performance between 
the shorter (e.g., Section 4.1.2) and longer time late values, detection probability drops by about 
30% for search missions lasting a single hour, but only shows an 8% decrease in performance 
after four hours of search. 

There is also a shift in the crossover point where myopic search begins to outperform linear 
search, as shown in Fig. 4.4(b). Compared to the previous analysis with t = 24hrs, now with 
a longer time late, this reversal occurs later at nearly four hours, highlighting the diminished 
advantage of using myopic search. The overall detection performance is also noticeably lower 
for this challenging scenario. 



Figure 4.4: Search mission results for T = 36 hours, W = 100m, and for (a) U = 2%, and (b) U = 0.2%. 
Linear search appears to perform better than myopic search no matter the search mission duration 
when the drifting objects move relatively faster. For slow-moving objects, only after four hours does 
myopic search begin to perform marginally better. 

4.1.5 Varying Initial Object Locations 

Following the preliminary analysis for the target starting at a point above the generalized source- 
sink pair, the target was made to start its drift from another point that is located below the 
generalized source-sink pair (see Fig. 4.5). The outcome of the simulations demonstrated that 
the trends described in the previous sections appear to hold for arbitrary starting locations of 
the drifting object. These corroborating results further reinforce the relationships between these 
parameters described previously, and highlight the merit of further, more extensive sensitivity 
analyses. 

4.2 Statistical Analysis on Potential Flow Model 

Using the Student’s T-test method, a hypothesis testing was carried out on the results. The null 
hypothesis, Hq, is defined as “both myopic search and linear search have the same detection 
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Figure 4.5: An example instance of an alternate initial location of the drifting object, such as below the 
source-sink pair. The functional dependencies of the probability of detection on the factors presented 
in this paper appear to hold for arbitrary initial locations. 

performance” and the alternative hypothesis, denoted as H\, is that the “detection performance 
between myopic search and linear search are distinctively different.” For the result to be deemed 
significant, we chose the ct-value to be 0.05 or 5%. 

For sweep width (IV) equal to 1000 m, we used the data analysis package in Microsoft Excel 
and performed a “T-test: Paired Two Sample for Means" on all the simulation results obtained 
from the various permutations of time late (t), search time (T) and drifter speed (U). The 
calculated p-value of 0.74 is significantly greater than a and thus we do not reject the null 
hypothesis (i.e., there is no evidence to infer that the alternative hypothesis is true). Next, we 
repeat the data analysis for W = 500 m. The calculated p-value was 0.0318, which is deemed 
significant as it lies between 0.01 and 0.05. Given that the p-value is less than a, there is strong 
evidence to infer that the alternative hypothesis is true and we reject the null hypothesis (i.e., 
the detection performance between myopic search and linear search are likely to be different). 
As expected, when IV is reduced further to 100 m, the calculated p-value was 0.0029, which is 
highly significant and provides overwhelming evidence to infer that the alternative hypothesis 
is true. 

The statistical analysis showed that, given a large W, the search performance between the two 
strategies will likely yield the same results. That is because the area created by the convex 
hull of the 500 samples is relatively small as compared to the area that can be covered by the 
searcher with a large TV of 1000 m. Intuitively, the ratio between the two areas starts to increase 
if IV is reduced by 50% to 500 m or 90% to 100 m and we are able to note a distinct difference 
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in detection probability results between linear search and myopic search. 

Fig. 4.6(a) shows the comparison of the mean detection probability and standard error for both 
myopic and linear search. For W = 1000 m, the mean and standard error obtained for myopic 
search is 99.19% and 0.29%, respectively. For linear search, the mean was 99.1% and a standard 
error of 0.15%. For W = 500 m, myopic search achieved a mean of 94.03% with a standard 
error of 1.15%, while linear search has a slightly higher mean of 95.43% and a standard error 
of 0.68%. For W = 100 m, the mean detection probability for myopic search was 38.8% with 
a standard error of 2.98%, while linear search continues to fare better with a higher mean of 
43.08% and a standard error of 2.33%. Notwithstanding, we note that there is a consistent 
sharp drop in detection probability when W is reduced from 500 m to 100 m, as shown in 
Fig. 4.6(b). This observation highlights the possibility of a non-linear relationship between 
detection probability and the linear decrease of sweep width. As such, additional values of W, 
spreading across 500 m to 100 m, needs to be investigated. Another interesting analysis is that 
the graphs in Fig. 4.6 showed that the naive linear search starts to outperform the myopic search 
for shorter sensing distances. 



Figure 4.6: Comparison of the mean probability of detection between myopic search and linear search 
for W of 1000 m, 500 m and 100 m 


4.3 Search Analysis Using Historical Data Set 

Historical data sets obtained from SWAFS provide the surface current speeds ( water u and 
water v variables) at hourly intervals for the months of March to August 2006. The hourly sur¬ 
face current data are spread out on a grid of 452 x 404 points covering the entire Persian Gulf 
region from 23.48°N, 47.49°E to 30.5°N, 57.51°E. For this thesis, we focused the simulation 
work and analysis primarily on the northern Persian Gulf region (28.5°N, 48°E to 30°N, 50°E) 
in order to draw references from the study done by Williams in 2007 on the High/Low wind 
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events during summer and winter. Similar to the potential flow model discussed in the earlier 
section, we chose two different start locations for the drifter as plotted in Fig. 4.7. Drifter start 
location 1 is about 100 km offshore from the nearest land mass on the west, while location 2 
is chosen to be about 30 km away from the northern land mass. The resultant separation of the 
two locations is approximately 88 km and the a priori understanding is that the effects of LSBC 
and tides should have a more pronounced effect at location 2 as compared to location 1, thereby 
generating stronger surface current speeds. 


Target starling points in North Arabian Gulf 



Figure 4.7: Drifter starting locations in north Persian Gulf. Location 1 is situated at 29°N, 49.2°E 
while location 2 is at 29.5°N, 48.5°E 

Using the 2D interpolation method in Matlab, the surface current speed and angle are calculated 
and the drifter position is recorded at intervals of five minutes. So, for t = 24 hours, there are a 
total of 288 data point calculations to show the drifter movement. This process is performed for 
drifters starting at location 1 and 2 as well as for the two dates chosen (i.e., 15 April and 15 July) 
in this analysis. The mean wind speed on 15 April is less than 5 m/s, whereas the mean wind 
speed on 15 July is higher than 5 m/s. Fig. 4.8 shows the drifter plots, without noise added, for 
T = 24 hours. While the quiver plot snapshots taken at the end of 24 hours for all four plots do 
point to the fact that surface current flows for both locations are generally in the same direction 
(i.e., on 15 April, the surface current speed is flowing towards northwest while on 15 July, the 
flow is to the east), the drifter plots showed that there is a distinct difference between the surface 
currents affecting the drifter movement at location 1 as compared to location 2, with location 2 
registering a higher mean target speed. 

The plot of the drifter positions found in Fig. 4.8 also portrayed the fluctuations in surface 
currents due to variations in tides and wind direction over a period of 24 hours. Taking an 
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Enlarged View of Target Location 1 
Mean Target speed 0.53 % of Searcher Speed, 


Enlarged View of Target Location 2 
Mean Target speed 0.85 % of Searcher Speed, 



(a) 


Enlarged View of Target Location 1 
Mean Target speed 0.53 % of Searcher Speed, 



(c) 



(b) 


Enlarged View of Target Location 2 
Mean Target speed 0.93 % of Searcher Speed, 



(d) 


Figure 4.8: Drifter plots, without noise, for T = 24 hours, on 15 April for (a) location 1, and (b) 
location 2, as well as on 15 July for (c) location 1 and (d) location 2. The magenta line represents 
the drifter movement and the * symbol is used to denote the final position of the drifter at the end 
of 24 hrs. 


hourly snapshot of the surface currents on 15 April, we found that at about T = 4 hrs, the 
surface current flow in the northern Persian Gulf started to move southeast along the major axis 
with an intensification of the flow at around T = 6 to 8 hrs before subsiding at T = 10 hrs and 
finally reversing the flow towards northwest at about T = 12 hrs (See Fig. 4.9 (a) to (f)). This 
flow pattern is repeated at semidiurnal periods with the same flow towards northwest at about 
T = 24 and 36 hrs (See Fig. 4.10 (a) and (b) respectively), suggesting that the dominant driver 
of the surface current speed on the 15 April is due to tidal currents. 

To illustrate the drifter speed, we determined the mean speed that the drifter was traveling and 
presented it as a percentage of searcher speed. For example, in Fig. 4.8(a) the mean target 
speed at location 1 on 15 April was 0.53% of searcher speed or in other words, with a searcher 
speed of 50 m/s, the target was traveling at an average speed of 0.265 m/s. That also means 
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Surface current flow of North Arabian Gulf starting from 15-Apr-2006 08:00:00 
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(e) (f) 

Figure 4.9: Surface current flow on 15 April, for (a) X = 4 hrs, (b) X = 6 hrs, (c) X = 8 hrs, (d) 
X = 10 hrs, (e) X = 11 hrs and (f) X = 12 hrs. Surface currents tend to flow towards the northwest 
at semidiurnal periods 


that the distance traveled by the drifter over a duration of 24 hours was 22,896 m. However, 
this value should not be interpreted as saying the drifter is now approximately 22 km away 
from where it first started. By inspecting the plots in Fig. 4.8, we are able to discern that 
the distance separating the start and final position of the drifter is actually much less than the 
distance traveled based on mean speed calculations. As such, if we had used the mean speed 
value and the mean direction of travel to plot the drifter position, it would result in a large error 
(as shown in Fig. 4.11) between the exact drifter position and the calculated mean position. 
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Surface current flow of North Arabian Gulf starting from 15 Apr 2006 08:00:0 
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Figure 4.10: Surface current flowing northwest on 15 April, for (a) z = 24 hrs and (b) z = 36 hrs 
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Figure 4.11: Illustration of the difference between using the mean speed (0.53% or 0.265 m/s) and 
direction (035°) to plot the drifter position in time as compared to using the interpolated speed and 
direction of the surface currents calculated at five-minute intervals. 


With the U speed being defined by the historical data set, the same factors and values that was 
used in the potential flow model was chosen for the simulation runs. One of the findings from 
the potential flow model was that there exists a possible non-linear drop in detection probability 
when W is reduced from 500 m to 100 m and greater resolution of W was needed. Assuming a 
linear drop in detection probability, we added additional values of W (i.e., 600 m, 400 m, 300 m 
and 200 m) that are spaced 100 m apart. Subsequent analysis showed that the steepest gradient 
occurred when W was reduced from 200 m and 100 m. As such, three additional values of W 
spaced 25 m apart were introduced into the simulation runs, resulting in a total of ten different 
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Table 4.2: List of factor values that are used in the simulation runs with drifter speed, U, being 
extracted from the historical data set provided by SWAFS. 


Factors 

Values Used in Simulation 

Time Late (t) 

24 and 36 hrs 

Search Time ( T) 

1, 2, 3 and 4 hrs 

Sweep width (W) 

100, 125, 150, 175, 200, 300, 400, 500, 600 and 1000 m 


values of W investigated with a total of 320,000 simulation runs conducted. Table 4.2 lists the 
different factor values that were investigated with the historical data set. 

As explained in Section 3.2.3, we introduced Weibull-distributed perturbations to the drifter 
paths at each time step to create random dispersions of the drifter cloud. The enclosing convex 
hull (i.e., formed by the dotted red line) shown in Fig. 4.12 thus correspond to the 3a expanding 
uncertainty area created by the 500 simulated particles. Thereafter, linear search and myopic 
search strategies (See Section 3.3.2) were performed on all 500 particles and the search results 
tabulated and expressed as a percentage of the total detection probability. Shown in Fig. 4.12 
are the search results for particle number 1 at varying sweep widths of 1,000 m, 500 m, 200 m 
and 100 m. In essence, the searcher is able to locate particle number 1 in all four illustrations 
but the time taken to conduct the search increases as W is reduced. If the search time, T, is 
limited to only the 1st hour, then the searcher is only able to locate the target if it uses a sensor 
with a sweep width of 1,000 m. Similarly, if T is increased to 2 hours, the minimum W setting 
must be set at 500 m in order to localize the target within the allocated search time. 

Another point to note and is also shown in Fig. 4.12(a) and (b) is that for W > 500 m, they 
will have the same Search Plan (i.e., [1, 3, 2, 4]). That is to say, the searcher will proceed to 
conduct its search operations in area 1 for the 1st hour, area 3 for the 2nd hour, area 2 for the 
3rd hour and finally, area 4 at the 4th hour. However, when W is reduced to less than or equal 
to 200 m, there is a shift in the order of the Search Plan with area 4 being assigned as the area 
to search for the 2nd hour, area 3 to be searched in the 3rd hour and area 2 to be searched in 
the 4th hour. This shift in search plan is a result of the “greedy” algorithm in myopic search 
technique as described in Section 3.1.2. With the reduction in sweep width, the searcher can 
no longer fully cover the entire area of uncertainty within the user-defined search duration of 1 
hour and thus needs to adopt a locally optimal plan that will maximize the chances of detecting 
the drifter while minimizing the mean time it takes to do it. 
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Weibull: Enlarged View of Spiral Search 
Mean Drifter speed 0.52916 % of Searcher Speed, Searcher Sweepwidth 1000 m 



Weibull: Enlarged View of Spiral Search 
Mean Drifter speed 0.52916 % of Searcher Speed, Searcher Sweepwidth 500 m 
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Weibull: Enlarged View of Spiral Search 
Mean Drifter speed 0.52916 % of Searcher Speed, Searcher Sweepwidth 200 m 


Weibull: Enlarged View of Spiral Search 
Mean Drifter speed 0.52916 % of Searcher Speed, Searcher Sweepwidth 100 rn 
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Figure 4.12: Search results for particle run number 1 starting at location 1 on 15 April for t = 24 
hrs with a variation of W ((a)1000 m, (b)500 m, (c)200 m and (d)100 m) being evaluated. As W 
reduces, the “greedy” algorithm recommends an alternative search plan that maximizes the detection 
probability while minimizing the search cost (i.e., time) 


Next, we present the results for 12 cases where the impact of the search and ocean model 
parameters on the probability of detection is evaluated. Cases 1 to 4 looked at the differences in 
detection probability between location 1 and 2 on 15 April 2006 when different values of time 
late (i.e., 24 and 36 hours) and sweep width (i.e., 500 and 100 m) are used. Similarly, cases 5 
to 8 looked at the results obtained on 15 July 2006 with the same permutations of time late and 
sweep width. Finally, cases 9 to 12 focused on the change in the detection probability curves 
for the entire sweep width range (i.e., 1,000 m to 100 m) as well as differnt values of time late 
for each location on 15 April. A summary of the different values used in each case is illustrated 
in Table 4.3. 
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Table 4.3: List of 12 cases and their associated parameters used in the simulation runs based on the 
historical data set. 


Case 

Date 

Time late (hrs) 

Sweep width (m) 

Remarks 

1 

15 Apr 

24 

500 

Comparing the difference in detection 

2 

15 Apr 

36 

100 

probability between location 1 and 2 

3 

15 Apr 

24 

500 


4 

15 Apr 

36 

100 


5 

15 Jul 

24 

500 


6 

15 Jul 

36 

100 


7 

15 Jul 

24 

500 


8 

15 Jul 

36 

100 


9 

15 Apr 

24 

1,000 to 100 

Detection probability curves 

10 

15 Apr 

36 

1,000 to 100 

at location 1 

11 

15 Apr 

24 

1,000 to 100 

Detection probability curves 

12 

15 Apr 

36 

1,000 to 100 

at location 2 


4.3.1 Case 1: 15 April, Time Late, t = 24 hrs; Sweep Width, W = 500 m 

Based on a time late of T = 24 hours, and a searcher with a sensor sweep width W = 500 meters, 
the simulation results for the case where the drifter was dropped at location 1 on 15 April is 
shown in Fig. 4.13(a), whereas Fig. 4.13(b) shows the results where the drifter was dropped at 
location 2. The mean speed at location 1 was 0.265 m/s (0.53% of searcher speed, V ) while the 
mean speed at location 2 was 0.425m/s (0.85% of V). The linear search approach continued to 
perform better than myopic search during the first search hour and appeared to plateau after two 
hours of search effort. This is because the drifting object had already moved out of the search 
zone within those two hours such that continued searching in subsequent areas was fruitless. 
We observed a crossover point where, although the linear search approach initially excelled, 
myopic search began to perform slightly better after search duration time of about two hours. 

While the drifter at location 2 was experiencing a mean speed that was 60% more than location 
1, the simulation results showed that it has relatively little impact when the searcher is using 
a sweep width of at least 500m. This fact again reemphasized the earlier point as proven in 
the potential flow model that, for small enough relative flow speeds, the expanding area does 
not increase too rapidly to be completely covered by a searcher with a sweep width of at least 
500m. 
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Figure 4.13: Probability of detections from simulation studies comparing linear search (red squares) 
and myopic search (blue diamonds). The search and ocean parameters are z — 24 hours, VF=500 m, 
and mean surface current speed or object speed ratio at location 1 (a) C/=0.265m/s (0.53% of V) 
and at location 2, (b) [/=0.425m/s (0.85% of V). For both cases, there is a crossover point at T — 
2hrs where myopic search performs slightly better than linear search 


4.3.2 Case 2: 15 April, Time Late, t = 24 hrs; Sweep Width, W = 100 m 

When the sensor sweep width was reduced from 500m to 100m, there was a sharp drop of about 
60% in the detection probabilities for linear search in both location 1 and 2. Myopic search 
also suffered from approximately 60% drop in location 1 and about 70% drop in detection 
probability in location 2. As illustrated in Fig. 4.14(a), with a slower mean surface current 
speed of 0.265m/s, there was a crossover seen only after search time, T, of four hours where 
myopic search began to outperform the linear search. However, when the mean surface speed 
was increased by 60% in location 2, there was no such crossover observed and linear search 
outperformed myopic search for all values of T (See Fig. 4.14(b)). 
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Figure 4.14: Search mission results for T = 24 hours, W = 100 m, and for mean surface current speed 
or object speed ratio at location 1 (a) [/=0.265m/s (0.53% of V) and at location 2, (b) [/=0.425m/s 
(0.85% of V). For shorter sensing distances and slower surface currents as in (a), the myopic search 
approach appears to only outperform the naive linear search at T = 4 hours 
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4.3.3 Case 3: 15 April, Time Late, t = 36 hrs; Sweep Width, W = 500 m 

As with the potential flow model, we next investigated the impact of increased time late of 
T = 36 hours (with W = 500 m), which nominally only results in a larger search area due 
to further expansion over time as compared with the above cases. The overall mean surface 
current speed increased by about 3% for location 1 to 0.2725 m/s (0.545% of V ) and about 
1.5% for location 2 to a mean speed of 0.4315 m/s (0.863% of V). 

In doing so, we observed that there was approximately 15% and 6% drop in detection probabil¬ 
ity for myopic search and linear search respectively, when we compared the results for location 
1 and 2 at T = 1 hour with the results in Section 4.3.1. However, this disadvantage was quickly 
regained, as both myopic search and linear search was able to produce similar search perfor¬ 
mances as in Section 4.3.1, when search time allocated to the mission was increased to two 
hours or more as illustrated in Fig. 4.15(a) and (b). What this means is that while the expand¬ 
ing search area may have grown bigger with T = 36 hours, it had twice the impact on myopic 
search than linear search. However, this degradation was only limited to T = 1 hour, as both 
search techniques were able to close the performance gap for missions extending past the first 
hour. In addition, the 58% difference in U between location 1 and 2 did not affect the overall 
performance of the mission between the two search techniques past the first hour. 



Figure 4.15: Search mission results for T = 36 hours, W = 500 m, and for mean surface current speed 
at location 1 (a) [/=0.2725m/s (0.545% of V) and at location 2, (b) [/=0.4315m/s (0.863% of V). 
For larger search area, myopic suffered more initially but managed to regain its performance past the 
initial hour. The 58% difference in U had minimal impact on search performance 


4.3.4 Case 4: 15 April, Time Late, t = 36 hrs; Sweep Width, W = 100 m 

As before, we examined what happened when the sensor sweep width was reduced to W = 100 
m, for extended time late T = 36 hrs. Relative to previous studies, we anticipated diminished 
search performance. As shown in Fig. 4.16 (a) and (b), for U = 0.2725m/s (0.545% of V) and 
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U = 0.4315m/s (0.863% of V ), respectively, we observed that there was about 17% drop in the 
highest probability of detection for myopic search when compared with the former case, even 
with an extended search mission duration of four hours. Linear search fared much better with 
only about 10% drop in the highest probability of detection when compared with the former 
case. Linear search also continued to outperform myopic search for all search times, as shown 
in Fig. 4.16 (a) and (b). 


Given that location 2 had a 58% increase in U, the overall detection performance was expected 
to be lower as compared to location 1. From Fig. 4.16, there was an approximately 10% differ¬ 
ence in search performance at the extended search duration of four hours for myopic search and 
about 6% for linear search when comparing the results between (a) and (b). 



Figure 4.16: Search mission results for T = 36 hours, W = 500 m, and for mean surface current 
speed at location 1 (a) C/=0.2725m/s (0.545% of V) and at location 2, (b) f/=0.4315m/s (0.863% 
of V). Linear search outperformed myopic search regardless of the search mission duration and mean 
U speed 


4.3.5 Case 5: 15 July, Time Late, t = 24 hrs; Sweep Width, W = 500 m 

Next, we looked at the surface current conditions on 15 July where the mean wind speed data 
points was above 5m/s. Based on a time late of t = 24 hours, and a searcher with a sensor 
sweep width W = 500 meters, the simulation results for the case where the drifter was dropped 
at location 1 and location 2 on 15 July is shown in Fig. 4.17(a) and (b) respectively. The mean 
speed at location 1 was 0.265 m/s (0.53% of searcher speed, V) while the mean speed at location 
2 was 0.463m/s (0.926% of V). That is to say, location 2 experienced a mean surface current 
speed that was 75% more than location 1. 

Linear search continued to perform better than myopic search during the first search hour and 
appeared to plateau after two hours of search effort because the drifting object had already 
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moved out of the search zone within those two hours. We observed a crossover point where 
myopic search began to perform slightly better after search duration time of about two hours. 
Similar to Section 4.3.1, while the drifter at location 2 was experiencing a mean speed that was 
75% more than location 1, the simulation results showed that it had relatively little impact for 
t = 24 hours and when the searcher was using a sweep width of at least 500m. The expanding 
area did not increase too rapidly and could easily be completely covered by a searcher with a 
sweep width of at least 500m. 



(a) 



(b) 


Figure 4.17: Probability of detections from simulation studies carried on two locations on 15 July. The 
search and ocean parameters were T = 24 hours, W=500 m, and mean surface current speed or object 
speed ratio at location 1 (a) f/=0.265m/s (0.53% of V) and at location 2, (b) [7=0.463m/s (0.926% 
of V). For both cases, there was a crossover point at T = 2hrs where myopic search performed 
slightly better than linear search 


4.3.6 Case 6: 15 July, Time Late, t = 24 hrs; Sweep Width, W = 100 m 

When the sensor sweep width was reduced from 500m to 100m, linear search experienced a 
drop of about 60% in detection probabilities for the first hour of search in both location 1 and 2. 
Myopic search, however, suffered from an average of about 65% drop in detection probabilities 
in location 1 and 2 during the first search hour. For both locations, there existed a crossover 
point after search time, T, of 3 hours where myopic search started to outperform linear search 
by about 10 to 12% at the extended search duration of four hours. Next, we compared the 
impact of the 75% increase in mean surface speed at location 2. We observed that there was a 
slight reduction of about 6% in detection probability, as shown in Fig. 4.18, for both myopic 
and linear search at the extended search duration of four hours. 

4.3.7 Case 7: 15 July, Time Late, % = 36 hrs; Sweep Width, W = 500 m 

We continued the same analysis and investigated the impact of increased time late of r = 36 
hours (with W - 500 m). The overall mean surface current speed decreased by about 4% for 
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Myopic 



(b) 


Figure 4.18: Search mission results for t = 24 hours, W = 100 m, and for mean surface current speed 
or object speed ratio at location 1 (a) U= 0.265m/s (0.53% of V) and at location 2, (b) £/=0.463m/s 
(0.926% of V). Myopic search approach appeared to only outperform the naive linear search at T = 
3 hours 


both locations 1 and 2. The mean surface current speed at location 1 is now 0.255m/s (0.51% 
of V) while location 2 is experiencing a mean speed of 0.443m/s (0.885% of V ). 

With a larger search area due to the increased time late, we observed that myopic search ex¬ 
perienced a drop in detection probability of about 16% and 18% at location 1 and location 2 
respectively, when we compared the results for location 1 and 2 at T = 1 hour with the results 
in Section 4.3.5. Similarly, linear search also suffered a drop in detection probability of about 
8% for both locations during the first hour of search. As such, the expanding search area with 
T = 36 hours has twice the impact on myopic search than linear search. 

When the search mission was extended to two hours or more as illustrated in Fig. 4.19, both 
myopic search and linear search were able to achieve comparable search performance with an 
average of about 98% in detection probability. In addition, the 74% difference in U between 
location 1 in Fig. 4.19(a) and location 2 in Fig. 4.19(b) only affected the overall performance 
of the two search strategies by a maximum reduction of 4% during the first search hour. The 
impact of the increased U was subsequently reduced to a maximum of about 1.6% when the 
mission was extended to four hours. 


4.3.8 Case 8: 15 July, Time Late, t = 36 hrs; Sweep Width, W = 100 m 

Next, with the sensor sweep width reduced to W = 100 m and for extended time late t = 36 hrs, 
we expected a diminished search performance. Fig. 4.20(a) and (b) show the simulation results 
for location 1, with U = 0.255m/s (0.51% of V) and location 2, with U = 0.443m/s (0.885% of 
V), respectively. When compared with the simulation results in Section 4.3.6, we observed that 
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Figure 4.19: Search mission results for T = 36 hours, W = 500 m, and for mean surface current speed 
at location 1 (a) U=0.255m/s (0.51% of V) and at location 2, (b) U=0.443m/s (0.885% of V). 
Myopic suffered more initially but managed to regain its performance past the initial hour. The 74% 
increase in U had minimal impact on the search performance 


there was about 13% drop in the highest probability of detection for myopic search for search 
mission duration of four hours. 

Linear search fared much better, with only about 2% to 5% drop in the highest probability 
of detection. Linear search also continued to outperform myopic search for search times less 
than or equal to three hours, as shown in Fig. 4.20(a) and (b). In addition, there existed a 
crossover point for both locations where, after a search time of four hours, myopic search began 
to perform slightly better than linear search. While location 2 had a 74% increase in U, the 
overall detection performance between Fig. 4.20(a) and (b) was about 5% difference in search 
performance at the extended search duration of four hours for both myopic search and linear 
search. 


4.3.9 Case 9: 15 April, Location l,Time Late, t = 24 hrs; Search Time, 
T = 1 to 4 hours 

Earlier simulations results from the potential flow model indicated the possibility of a non-linear 
relationship between sweep width and detection probability. However, as there were only three 
settings of sweep width (i.e., 100, 500 and 1,000 m) used in the potential flow model, it was 
difficult to discern if the relationship between sweep width and detection probability was linear 
or otherwise. As such, we added a total of seven more sweep width settings into our simulation 
runs as shown in Table 4.2 and the simulation results showed that the detection probability for 
both linear and myopic search strategies decreased non-linearly with sweep width, with the 
sharpest drop in the detection probability of about 30% to 40% between 200 m and 100 m. 
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Myopic 



(b) 


Figure 4.20: Search mission results for t = 36 hours, W = 500 m, and for mean surface current speed 
at location 1 (a) £/=0.255m/s (0.51% of V ) and at location 2, (b) C/=0.443m/s (0.885% of V). 
Linear search outperformed myopic search for search mission duration less than or equal to 3 hours. 
Myopic search only managed to perform slightly better when the search duration was extended to 
four hours. The 74% increase in U caused a 5% drop in detection probability 


Further analysis in Section 4.4 will provide an explanation of the cause of this non-linearity 
relationship. 

For now, we focused the analysis on the positive impact of extending mission search time, 
T, as well as the response between myopic and linear search on increasing T. As shown in 
Fig. 4.21(a), for T = 1 hour, we observed that linear search performed better for all ranges 
of sweep width. When T was increased to two hours, myopic search was able to achieve 
comparable search performance and performed slightly better than linear search when the sweep 
width is more than or equal to 400 m. There existed a crossover point at W between 400 and 
300 m as shown in Fig. 4.21(b), where linear search performed better than myopic search. 
Similar detection probability curves were observed for T = 3 hours and 4 hours as illustrated 
in Fig. 4.21(c) and (d). We observed that the crossover point where linear search started to 
perform better was progressively shifted towards a sweep width of 300 m as the search duration 
was increased further. By extending the mission search time to four hours, we also observed an 
approximately 10% increase in detection probability when W was set at 100 m. 


4.3.10 Case 10: 15 April, Location l,Time Late, t = 36 hrs; Search Time, 
T = 1 to 4 hours 

When T was increased to 36 hours, we noted from our earlier analysis that the increase in 
search area uncertainty would lead to a reduction of detection probability. The drop in detection 
probability was more pronounced for myopic search when the search duration was set at 1 
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Figure 4.21: Search mission results at location 1 for T = 24 hours, W ranging from 100 m to 1,000 
m, and for varying search time with (a) T = 1 hour, (b) T — 2 hours, (c) T — 3 hours and (d) 
T — 4 hours. Linear search outperformed myopic search for all sweep width settings at search mission 
duration equal to 1 hour with myopic search performing better at W > 400 m for extended durations. 
The crossover point where linear search performed better was observed to progressively shift towards 
W of 300 m when T increased 


hour as shown in Fig. 4.22(a). When T was increased to two hours, linear search continued 
to perform better than myopic search at W between 400 and 300 m as shown in Fig. 4.22(b). 
Similar detection probability curves were observed for T = 3 hours and 4 hours as illustrated 
in Fig. 4.22(c) and (d). With increased mission search time of four hours, we observed that 
myopic search performance started to resemble that of linear search. 


4.3.11 Case 11: 15 April, Location 2, Time Late, t = 24 hrs; Search Time, 
T = 1 to 4 hours 

Next, we looked at location 2 where the mean surface current speed affecting the drifter was 
60% more than location 1. For T = 1 hour, linear search continued to perform better for all 
ranges of sweep width as shown in Fig. 4.23(a). When T was increased to two hours, the 
crossover point where linear search started to perform better than myopic search shifted to 
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Figure 4.22: Search mission results at location 1 for t = 36 hours, W ranging from 100 m to 1,000 
m, and for varying search time with (a) T = 1 hour, (b) T = 2 hours, (c) T = 3 hours and (d) T = 4 
hours. Myopic search performed considerably worse than linear search at T = 1 hour but started to 
resemble linear search performance as T was increased to 4 hours 


around W of 500 m. As mission duration increased further, the crossover point was shifted down 
to a W that was between 400 and 300 m as shown in Fig. 4.23(b), (c) and (d). In addition, we 
observed that when T was increased to three hours, myopic search had about 10% improvement 
in detection probability at W of 300 m but minimal improvement for W at 200 m and below. On 
the contrary, the improvement of detection probability for linear search was more gradual and 
reaching its plateau towards T = 4 hours. 

4.3.12 Case 12: 15 April, Location 2,Time Late, t = 36 hrs; Search Time, 
T = 1 to 4 hours 

Similar analysis can be seen from Fig. 4.24(a) as we note that myopic search continued to per¬ 
form much worse than linear search at T = 1 hour but quickly regained the search performance 
when T was increased to two hours. Crossover point was again seen at W between 400 and 
300 m as illustrated in Fig. 4.24(b). Another point to note is that the improvement in myopic 
search for W at 300 m and below was seen to improve, with the highest improvement of about 
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Figure 4.23: Search mission results at location 2 for T = 24 hours, W ranging from 100 m to 1,000 
m, and for varying search time with (a) T = 1 hour, (b) T = 2 hours, (c) T = 3 hours and (d) 7 = 4 
hours. Linear search performed better at T = 2 hours at W of about 500 m. This crossover point 
shifted down to a W of between 400 and 300 m as T increased to 3 hours 


10% registered for W at 200 m when comparing between Fig. 4.24(b) and (c). As T was further 
increased to four hours, the improvement of detection probability for W less than or equal to 
200 m was marginal, indicating that both linear search and myopic search were starting to reach 
their respective plateaus in performance. 


4.4 Effects of Weibull Probability Distribution Function 

Throughout the analysis of both the potential flow model and historical data sets, the noise data 
generated from the Weibull PDF was fixed at a scale factor A = 1.0 and a shape parameter, 
y = 1.5 for each coordinate as described in Section 3.2.3. In addition, the analysis carried out 
in Sections 4.2 and 4.3 points to a nonlinear relationship between the detection probability and 
the reduction of sweep width from 300 m down to 100 m. The hypothesis was, therefore, that 
if the 500 particles are spread farther apart from their mean value (i.e., by increasing the scale 
factor), there would be a proportionate reduction in detection probability. As such, we generated 
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Figure 4.24: Search mission results at location 2 for t = 36 hours, W ranging from 100 m to 1,000 
m, and for varying search time with (a) T = 1 hour, (b) T = 2 hours, (c) T = 3 hours and (d) 7 = 4 
hours. Linear search continued to dominate at T = 1 hour. Same crossover point was observed at W 
between 400 and 300 m for T greater than 1 hour 


a new set of noise data by increasing the scale factor A by 10% to 1.1 while keeping the shape 
parameter, y, unchanged. Next, we used the new noise data and ran the simulations for the 
drifter starting at location 1 with a sweep width ranging from 300 m to 100 m and a time late, 
T = 24 hours. The results obtained validated the relationship between particle spreading and 
detection probability. As illustrated in Fig. 4.25, a 10% increase in A was able to have an effect 
of reducing the overall detection probability by an average of 22.5% between the two data sets. 

Next, we used the “Trendline” function in MS Excel to determine if the plots show a linear 
or non-linear relationship. For search time T equal to 1 hour (See Fig. 4.26(a)), the detection 
probability curves for sweep width from 300 m to 100 m pointed to a linear relationship with 
a high R 2 value of 0.936 for the original myopic data and 0.983 for the myopic data with A 
increased to 1.1. Increasing T to two or three hours also indicated a linear relationship between 
sweep width and detection probability with a high R 2 value of more than 0.92. However, the 
linearity relationship was not evident when T was increased to 4 hours as shown in Fig. 4.26(b). 
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Figure 4.25: Comparison of detection probability between the original myopic data at location 1 with 
T = 24 hours and when the scale factor of the Weibull PDF was increased by 10% 


A quadratic function was used instead to produce a high R 2 value of 0.9821 for the original 
myopic data and 0.965 for the myopic data with X increased to 1.1. 

If we refer back to Fig. 4.21(a), we note the nonlinear relationship for T = 1 hour started at 
a sweep width of 400 m instead of 300 m. So, if we had used the results starting from a 
sweep width of 400 m, we would be able to obtain a quadratic function to explain the nonlinear 
drop in detection probability. For T = 4 hours, using the data set from sweep width of 300 
m onwards was sufficient to show the nonlinear drop in detection probability. As such, the 
nonlinear reduction in detection probability seen at lower sweep widths can be attributed to a 
confluence of factors used in the simulation, namely, the Weibull PDF parameters, the sensor 
sweep width as well as the amount of search duration available. 


4.5 Effects of High/Low and High/High Wind/Current Events 

References are made to the work carried out by Williams in 2007 in analyzing the High/Low 
Wind/Currents events in the Persian Gulf. For 15 April, it was characterized as a High Wind, 
Low Current event. In other words, the Northern Persian Gulf region was experiencing low 
surface currents but relatively strong offshore and alongshore winds on 15 April. This obser¬ 
vation was consistent with the mean surface current speed obtained from the SWAFS data set. 
For a time late, t = 36 hours, location 1 registered a mean surface current speed of 0.2725 
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Figure 4.26: Comparison of increasing the scale factor of the Weibull PDF at location 1 by 10% for 
(a) T = 1 hour and (b) T = 4 hours. Blue diamonds denote the original myopic search results with 
A = 1.0 while red squares denote the myopic search results with A = 1.1 


m/s based while location 2 registered a mean surface current speed of 0.4315 m/s. Given that 
location 2 was relatively nearer to land mass, it was more susceptible to strong offshore winds 
as compared to location 1. As shown in Subsection 4.3.3 and 4.3.4, the 58% difference in U 
had minimal impact on search performance when large sensor sweep width were used at the 
two locations. However, if the sweep width was reduced to 100 m, myopic search suffered an 
additional 17% drop in search performance due to the increased U. On the other hand, linear 
search only encountered a 10% drop in search performance. 

Similarly, 15 July was characterized as a High Wind and High Current event. The combine 
effects of high surface currents and high offshore winds created a mean surface current speed at 
location 2 (See subsection 4.3.7) that was 74% stronger than location 1. We observed the same 
minimal impact to search performance when the sweep width was large. However, when the 
sweep width was reduced to 100 m as described in subsection 4.3.8, the difference in search per¬ 
formance between both locations was only a reduction of about 5% in detection probability for 
both myopic and linear search. Thus, we noted that by using a relatively large sweep width, we 
were able to nullify the fluctuations in mean surface current speeds. But when the sweep width 
was reduced to 100 m, the effects of the wind and surface currents were not straight forward. 
Additional research work (e.g., run simulations for multiple start locations and different time 
periods) will be required to investigate and understand how the search performance is being 
affected by the High/Low Wind/Current events and the use of small sensor sweep widths. 
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Table 4.4: Calculated p-v alue result for the 10 different sweep widths used in the simulation runs. For 
each sweep width, there are 32 pairs of mean detection probability comparing between myopic search 
and linear search. Each mean detection probability result was obtained by performing 500 simulation 
runs. Blank entries denote p-values of < 0.00001, denoting statistically statistical difference between 
myopic and linear search methods. 


W( m) 

1000 

600 

500 

400 

300 

200 

175 

150 

125 

100 

p-value 

0.337 

0.089 

0.071 

0.062 

0.003 

- 

- 

- 

- 

- 


4.6 Statistical Analysis on Historical Data Set 

Similar to Section 4.2, we use the Student’s T-test method to carry out a hypothesis testing on the 
results obtained from the historical data set. The null hypothesis, Hq, remains unchanged and 
is defined as “both myopic search and linear search have the same detection performance.” The 
alternative hypothesis, denoted as H\, is defined as “the detection performance between myopic 
search and linear search are distinctively different.” For the result to be deemed significant, 
we chose the a-value to be 0.05 or 5%. The p-value result for all 10 sweep width values are 
calculated and listed in Table 4.4. We observed that for W > 400 m, the p-value is greater 
than a and thus we do not reject the null hypothesis (i.e., there is no evidence to infer that the 
alternative hypothesis is true). For W < 300 m, the calculated p-value is very much less than 
a. Given that the p-value is highly significant, there is overwhelming evidence to infer that the 
alternative hypothesis is true and we reject the null hypothesis (i.e., the detection performance 
between myopic search and linear search are likely to be different). The statistical analysis 
showed that as W gets larger, the search performance between the two strategies will likely 
yield the same results. This conclusion echoes the analysis carried out in Section 4.2. 

The mean detection probabilities for the 320,000 simulation runs performed on the historical 
data set are presented in Fig. 4.27. The plots reaffirmed the nonlinear decline in detection 
probability as the sweep width is reduced linearly. The nonlinearity starts at about W of 400 m 
to 300 m and, as described in Section 4.4, is closely related to a number of contributing factors. 
For example, if we increase the scale factor used in the Weibull PDF, it increases the spread of 
the 500 particles from their mean and thus reduces the overall detection probability. Given the 
same sweep width values and the amount of search time available, the nonlinearity will start at 
higher sweep width ranges. 

Next, we plot out the mean detection probabilities and their corresponding standard error bars 
as shown in Fig. 4.28. The immediate take away is that linear search has a better mean search 
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Figure 4.27: Total mean probability of detection between myopic search and linear search. The 
nonlinearity drop in detection probability for both search strategies was observed to start from a W 
of about 400 m to 300 m 


performance than myopic search for all values of W. For W > 400 m, there are overlaps 
between the two standard errors obtained for both linear and myopic search with myopic search 
registering an average standard error of 1.36% as compared to an average standard error of 
0.54% for linear search. For W < 300 m, there are no overlaps between the two standard errors 
even when myopic search and linear search registered a higher mean standard error of 3.15% 
and 1.99%, respectively. 
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Figure 4.28: Comparison of the total mean probability of detection with error bars for myopic search 
and linear search. The error bars overlap at W > 400 m. Even with a higher standard error for W < 
300 m, there are no overlaps between the two search strategies 
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CHAPTER 5: 

Conclusions and Future Work 


While today’s sea mines are more technologically advanced and stealthy, they are also more 
costly and require an in-depth knowledge of the various influence parameters and programming 
to allow them to strike the intended target at the opportune time. However, a rogue nation that 
has been sanctioned or one with little resources to bear against a superior navy is likely to lever¬ 
age cheap and low technology contact mines [35] that, when deployed through indiscriminate 
mine laying, will be able to quickly deter or deny enemy access to its littoral environments. 

The mine warfare mission has never been more important, with at least 30 countries having the 
ability to manufacture naval mines and about 20 of those exporting the mines they produce. All 
in all, more than 50 nations have attained a certain degree of mining capability and given their 
flexibility and cost-effectiveness, mines make attractive weapons in asymmetric warfare. The 
cost of producing and laying a mine is usually anywhere from 0.5% to 10% of the cost it takes 
to remove it. In addition, it can take up to 200 times as long to clear a minefield as compared 
to the time it takes to lay it. In summary, countering terrorist mining is very difficult and time 
consuming. However, combining the right tools (Ocean and Search models) and adding a little 
bit of innovation to it, a tough problem may be broken down to a much simpler problem that is 
much easier to solve. 


5.1 Simulation Model 

The search and localization of a drifting object have significant operational relevance in numer¬ 
ous contexts and the motivation for this topic was derived from a BAA released by ONR in 
2011 looking for innovative solutions in the localization and tracking of surface or near-surface 
drifting mines. Given the prevalence of unmanned systems as well as the advancement of 
computational resources and analytic methodologies, the stage is set for the development of in¬ 
tegrated models incorporating both oceanographic and search-theoretic methods. The presented 
work shows that a single unmanned system equipped with a capable detector can effectively re¬ 
spond to a time-delayed cue for an object under the influence of ocean surface currents. Monte 
Carlo simulations were carried out on a simplified ocean current model as well as historic ocean 
currents data obtained from SWAFS. The stochasticity of real-world ocean environments was 
addressed by introducing perturbations according to Weibull noises at each time interval, which 
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resulted in an expanding area enclosing the likely locations of the drifting object. Search of this 
simulated expanding area was conducted by using either a naive linear spiral-out search pattern 
or an optimized myopic allocation of search effort, and the probability of detection as a function 
of relevant parameters was investigated. 

5.2 Findings 

Various factors can affect the search outcome, and this study has shown that the speed of the 
surface current (relative to the search speed) and the sensor sweep width of the searcher have 
dominant influence on the probability of detection of the object. When the sweep width is large 
(i.e., greater than or equal to 400m), linear search and myopic search exhibit similar probability 
distribution curves. As the sweep width is reduced linearly down to 100m, we observed a 
nonlinear decrease in detection probability, which is attributed to the scale and shape factors 
used in the Weibull Noise. The use of large sweep width settings has also shown to be able to 
nullify the drop in detection probability brought on by an increase in surface current speeds. 
For a 58% to 74% increase in surface current speed, there was minimal change to the detection 
probability when the sensor sweep width used was large. However, when the sweep width 
was reduced to 100m, the drop in detection probability can range from 10% to 17% for a 58% 
increase in surface current speed to as little as 5% for a 74% increase in surface current speed. 
Further simulation runs will be required to fully map out the relationship between small sweep 
width and the change in surface current speeds. In essence, the presented work highlights 
the fact that for most cases, straightforward linear search usually performs better than myopic 
search in both detection probability and search mission time required to successfully find the 
drifting object. 

5.3 Areas of Future Research 

This study has clearly demonstrated the value proposition that ocean models and search-theoretic 
methods can bring in countering surface or near-surface drifting mine threat. While the naive 
linear search performed better than the optimized myopic search effort in this study, there are 
several other more advanced search algorithms (e.g., multiscale adaptive search [29]) that have 
been developed for both single and multiple searchers scenario. In addition, both ocean models 
and search Tactical Decision Aids are usually two separate and distinct models and additional 
efforts will be required to encapsulate these integrated models into a user interface for decision 
support software. This thesis focused on two-dimensional flow models to address the surface 
or near-surface mine threat. Future extensions that examine full three-dimensional flow mod- 
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els that include a myriad of surface and bottom mines would also enhance the viability of the 
presented approach. 
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